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A VASICEK-TYPE SHORT RATE MODEL WITH 
MEMORY EFFECT 

AKIHIKO INOUE, SHINGO MORIUCHI, AND YUSUKE NAKAMURA 


Abstract. We introduce a Vasicek-type short rate model which has 
two additional parameters representing memory effect. This model 
presents better results in yield curve fitting than the classical Vasicek 
model. We derive closed-form expressions for the prices of bonds and 
bond options. Though the model is non-Markov, there exists an as¬ 
sociated Markov process which allows one to apply usual numerical 
methods to the model. We derive analogs of an affine term struc¬ 
ture and term structure equations for the model, and, using them, we 
present a numerical method to evaluate contingent claims. 


1. INTRODUCTION 

The Vasicek model introduced by [12] is a classical short rate model. It 
is defined by a stochastic differential equation (SDE) of the form 

(1.1) dr{t) = {a — br{t)}dt + adW*(t) (t > 0) 

describing the short rate process {r(f)}, where a, b and a are positive 
constants and {IV*(f)} is a Brownian motion under the equivalent mar¬ 
tingale measure. In this paper, we introduce a Vasicek-type short rate 
model with memory effect which has some good properties. 

To define the Vasicek-type model, denoted by A"!, we replace the Brow¬ 
nian motion {W*(t)} in (II.Ilf by a stochastic process {Z(t)} which was 
introduced and investigated by mmm (see Definition I27TI below). The 
process {Z(t)} is a Gaussian, stationary increment process which has pa¬ 
rameters p, q representing memory effect. The process {Z(t)} is also an 
Ito process unlike, e.g., fractional Brownian motion. This implies that we 
can apply the standard stochastic calculus, that is, the Ito calculus, to the 
model M.. 

For the Vasicek-type model A4 and 0 < t < T, we derive a closed-form 
expression for the time t price P(t, T ) of a zero-coupon bond with maturity 
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T, which we also call a T-bond (see Theorem 13.II below, the proof of which 
is given in Section [I]) . From the result and the forward measure method 
(cf. nu and [H Chapter 7]), we also obtain a closed-form expression for 
the prices of European bond options (see Proposition 13.41 below!. 

The closed-form expression for P(t, T) also gives that for the yield 
Y(t,T) defined by 

(1.2) Y(t, T) = - 1 ° g ; P ( t ’ t r) (t<T) 

or P(t,T ) = exp[—(T — t)Y(t,T)] (see Theorem 13.21 below). As in other 
short rate models, we can estimate the parameters of the Vasicek-type 
model A4 by fitting the yield curve T i—>• Y (0, T ) to actual yield data. We 
will see that the model A4 shows better fitting results than the classical 
Vasicek model, thanks to the additional parameters p and q (see Section 

ED- 

The model A4 has memory effect. In other words, it is a non-Markov 
model. In general, non-Markov interest rate models have the drawback 
that numerical computations become difficult in them. It turns out that 
the model A4 is free from this drawback, thanks to the existence of an 
associated two-dimensional Markov process, which is obtained by coupling 
the short rate process (r(t)} with another process (see (15.1|) below). This 
associated Markov process allows one to apply usual numerical methods. 
In fact, using the Markov process, we derive analogs of an affine term 
structure and term structure equations (cf. [5) E] and [7, Chapter 5]) for 
the model A4, and present a numerical method to evaluate European 
contingent claims based on them (see Section [5} • 


2. The model 


Let {P t }t>o be the Q-augmentation of the filtration generated by a one- 
dimensional standard Brownian motion {W*(£)}t>o defined on a complete 
probability space (fl, P,Q). 


Definition 2.1. For real numbers p, q such that 

(2.1) 0 < q < oo, — q < p < oo, 

we define a process {Z(t)} t >o by 

(2.2) Z(t) = W*(t) ~ J {J pe- {p+q){s - T) l(T)dW*(r)^ ds 

where the positive deterministic function l is defined by 

2 qp 


(t > 0), 


l(r) = 1 


(p + 2q) 2 e 2 i T — p 2 


(r>0). 
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As is eaily seen, the process {Z(t)} is a continuous, Gaussian, (T)}- 
adapted, Ito process. It may not seem so, but {Z(t)} is also a stationary 
increment process. In fact, if {LL(f )}*<=« is another Brownian motion on 
another probability space (fhJ 7 , Q), then, for the Gaussian, stationary 
increment process {X(£)} t > 0 defined by 

(2.3) X(t) = W(t ) ~ f l f pe~ 

J 0 L J —oo 

the law of {Z(t)} t > 0 under Q is equal to that of {X(f)} t > 0 under Q. See [2, 
Theorem 5.2 and Example 5.3] and [TQ1 Section 2]); (12.3p corresponds to 
(28) in |T] or (1.1) in [TO] - If p = 0, then {Z(t)} reduces to the Brownian 
motion {W*(t)}. 

For a, 6, a G (0, oo), we consider the following Vasicek-type SDE: 

(2.4) dr(t) = {a — br(t)}dt + adZ(t), t > 0, r(0) G [0, oo). 

From 

d{e bs r(s )} = e bs dr(s ) + be bs r(s)ds = ae bs ds + cre bs dZ(s), 


^+i)( s ~^dW(r)\ds (t > 0), 


we see that the unique strong solution to (I2.4ji satisfies 

r(r) = e-h^rit) + ^(1 - e" 6 ^) + a [ e~ b{r - s) dZ(s) 
(2.5) W b K it 


(0 < t < r). 


We consider a Vasicek-type short rate model A4 in which the short 
rate process {r(t )} t >o follows (I2.4ji or (12.51) . We define the money-market 
account process {B(t)} t >o by 

B{t) ■.= eti r{s)ds (t > 0). 


Let T e (0, oo). For 0 < t < T, let P(t, T ) be the time t price of a T-bond. 
We regard Q as an equivalent martingale measure of A4 in the sense that 
the discounted price process 

(2.6) P( t ,T):=^p (0 <t<T) 

becomes a {Q, T)}-mart ingale. Then we have 


(2.7) 


P(t,T ) = E q 


e ~St r ^ s ) ds 



(0 <t< T ). 


For generalities of short rate models, one can consult, e.g., [7J Chapter 5]. 


Remark 2.2. If p = 0, then Z{t ) = W*(t), whence the Vasicek-type 
model M above reduces to the classical Vasicek model. 
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3. Prices of bonds and bond options 


Though the short rate process {r(t)} defined by (12.41) is not a Markov 
process, the model A4 admits a closed-form expression for P(t,T). To 
state this result, we introduce the following function: 

m{t) := / {l — j ds. 


We can write m(t) more explicitly as 


P 


+ 


bpe 


~(p+q)t 


pe 


—bt 


m(t) = < 


p + q (p + q - b)(p + q) p + q - b 

=-(p+<?)* 


(p + q-b^ 0), 


P 


pe 


— pte 


—bt 


(p + q- b = 0). 


p + q p + q 

Here is a closed-form expression for P(t,T). 

Theorem 3.1. For P(t,T) in (12.71) with (12.41) and (12.2D . we have 
P{t , T) = exp {-A(t, T ) - C(t, T)r(t) + U(t , T)} (0 < t < T), 

where 

1 _ e -b(T-t) 


C(t,T) := 


Mt,T) 


r T-t 


|m(s) + e bs — l} 2 ds 

a 2 qm 2 iT — t) 
b 2 {{p + 2 q) 2 e 2qt — p 2 } ’ 


U(t, T) := 


crip + 2 q) 2 e qt m(T — t ) 


3 (JS _ 


P 




i 


dZ(s). 


6{(p + 2q) 2 e 2 i t — p 2 } J o V p + 2g / 

We prove Theorem 13.11 in Section [H 

Recall Y(t,T) from (11.21) . From Theorem 13.11 we immediately obtain 
the next theorem. 

Theorem 3.2. For Y(t,T) in (11.21) with P(t,T) as in Theorem 1,7.71 we 
have 


ni , T) = 4M) + £M) r(() _^.D 


T — t 


T-t 


T-t 


(0 < t < T ). 


In particular, 

(3,1) = ^ + (T> 0). 

In Section [6j we use (13.ip to estimate the parameters of the model A4 
from actual yield data. 

Recall P(t,T) from (12.61) . 
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Proposition 3.3. The process {Pit, P)}o<t<T satisfies the SDE 

(3.2) dP(t, T ) = v(t, T)P(t, T)dW*(t) (0 < t < T), 
where the deterministic function v(t, T) is defined by 

(3.3) v(t,T) := ^{e" 6(T_t) -l + Z(t)m(r-t)} (0 <t<T). 

Proof. For fixed T > 0, we put 

V{t) := -Aft, T ) - C(t, T)r(t ) + U{t , T) — f r{u)du (0 <t<T), 

Jo 

so that P(t,T) = exp{V(t)}. Then, from ( 12 .4p and ( 12 . 21 ) . we have 


dV{t) = ${t)dt 

, „ f_^/, , (P + 2g) 2 e gt m(T - t) 

\ J 6((p + 2g) 2 e 2 ^-p 2 ) 

= P(t)dt + v(t,T)dW*(t) 



V 

p + 2q 



dW*(t ) 


for some continuous {Jy}-adapted process {fi{t)}o<t<T- This and Ito’s 
formula yield 

dP{t, T) = v(t, T)P(t, T)dW*(t) + 7 (t)dt (0 < t < T) 


for another continuous {T) {-adapted process { 7 (t)}o<t<T- However, since 
{P(t, r)}o<t<T is a martingale, we see that 7 (t) = 0. Thus the proposition 
follows. □ 


Let 7 Tc( 0 ) be the time 0 price of a European call option on a T-bond, 
with maturity S £ (0, T) and exercise price K £ (0, 00 ). Then 7 rc( 0 ) is 
given by 


(3.4) 


vr c (0) = E Q 


e ~ Jo r ( s ) ds (P(S,T) — K) + 


) 


where (x) + := max(a:,0). We write <f> for the cumulative distribution 
function of a standard normal distribution: 


$(x) := / j_ e <yl ^ 2 ' )y2 dy. 

K ! 7-oo VrU 


Proposition 3.4. For 7 r<y( 0 ) in (13.4|) with {r(£)} and {P{t,T)} as in 
Theorem, IT 11 we have 


7T C (0) = P(0, T)$(d + ) - KP( 0, S)$(d_), 
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where 


d i := 




log 


mr) 

KP{ 0, 5) 


+ ^\S) 


fS 

d _ := d+ - v/S%S), S 2 (5) := / v\ T (t)dt, 

Jo 


v S ,t(J) ■— ^ {e 6(T - e 6(5 4) + l(t)(m(T - t) - m(S - t))} . 


Proof. Notice that vs,t(1) = v(t, T ) — v(t, S ). Since v(t, T ) is a determin¬ 
istic function, the proposition follows immediately from Proposition 13.31 
and [3 Proposition 7.2], □ 


We use Proposition l3.4l to illustrate the efficiency of a numerical method 
in Section [5j 


Remark 3.5. The Vasicek model (11.11) is one of the affine term-structure 
models which are extensively used in finance. The CIR model is another 
example. See mum- See also 111 for related SDEs and diffusion 
processes. It is an open problem to obtain results similar to those in this 
paper for the CIR-type and other models that are obtained by replacing the 
Brownian motion by the process {Z(t)}. 


4. Proof of Theorem 13.11 


From (12.51) . 

rT 


r{r)dr = a 


e~ b{T - s) dZ(s) \ dr 


+ C(t 1 T)r(t) + -{T-t-C(t,T)}. 


By the classical Fubini theorem as well as that for stochastic integrals (cf. 
[3 Section 6.5]), the first term on the right-hand side is equal to 

r T 1- e -b(T-s) roo 


a 


-dZ(s ) = -a / f(s)dZ(s), 


where 


r e -b{T-s) _ x j 

/(*) : = --r- 


Combining these with (12.71) . we obtain 


(4.1) 


P(f T) = exp —- (T — t — C(t, T)) - C(t, T)r 


x 


E q 


exp < a / f(s)dZ(s ) 


T t 
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We write P\o,t] for the orthogonal projection from L 2 (f2, T,Q) onto the 
closed subspace spanned by (Z(r) : r G [0, £]}, and PW for its orthog¬ 
onal complement: P^X := X — P[o,t]X for X G L 2 (fi, X, Q). From [2j 
Theorem 5.2 and Example 5.3] and [TDJ Section 2], we have, for t > 0, 

^ +29) Ce±*£2j^ (t) U. 


W*(t) = Z(t) + / , , ^ 

7o Uo ' (p + 2q)e qT — p 2 e~ qT 

This and (12.21) imply 

(4.2) £r(Z(s) :«G[0,t]) = a(W*(s) : s G [0, t]) (t > 0). 

Since the process {Z(s)} is Gaussian, we see from (14.21) that T t and 
P[q t ] f t °° f(s)dZ(s) are independent. Hence, 


E q 


(4.3) 


exp < a / f(s)dZ(s ) 


P 


= exp aP[ 0 ,t] / f(s)dZ(s ) 


a 

x ex P ^ y 


P[o,t] / /(a)dZ(a) 


where ||X|| := P^x 2 ] 1 / 2 for X G L 2 (H,P,Q). 

First, we calculate P[ 0jt ] / t °° f(s)dZ(s). From [21 Theorem 4.7 and Ex¬ 
ample 5.3], 


P [0 , t] / /(a)dZ(a) = 


'0 wo 


s(s, t; 0, t)f(t + r)dr j dZ(s), 


where 


pe 


-gt 


0(f) : = , 9 , 

p + 2q 

b(s, r) := -pe- qs e~ (p+q)r , 

g(s , r; 0, t) := b(t - s, r)- -+ H 5 , r )‘ ^ 


1 — 0 2 (t) ' ' ' 1 — (j) 2 (t) 

However, since I[t.T] (t + r) = /[ 0i T_t](r), we have 

POO 

/ 0(s,T;O,t)/(t + T)dT 

J 0 

p + cf){t)e~ qs j 


rT-t 


H 1 - 0 2 O)} 

e - qt m(T — t ) 

K 1 - 0 2 W} 


e 9S - 


e -(p+q)r | e -6(T-t-r) _ 1 j rfr 

) 


p + 2g 


Thus 

(4.4) 


P[ 0 , t ] f°° f(s)dZ(s) = --U(t,T). 

Jt ° 
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Next, we calculate \\P^ t ] / t °° f (s)dZ (s)\\ 2 . By (2J Theorem 4.12 and 
Example 5.3], 


P[it] / f(s)dZ(a) 


oo 


n=0 ' 


d 2 n (s,f)ds, 


where 


do(sj) 


and 


c{u)f{t + s + u)du + f(t + s) 


(s > 0) 


(4.5) 


d n {s,f):=—J c('«) | J b n [t+ u+s,T)f(t+ T)drj du 

roo 

+ / b n (t + s, r)f(t + r)dr (n e N, s > 0) 


with 

c(s) := pe _(p+9)s J (0 ,oo)(s), 

£> n (s, r ) ;= -^-\ t )pe~ qs eP p+q)T (n e N, s, r > 0). 
Since I[ t ,T](t + s + u) = I[ 0tT -t- s ]{u)I[o, T -t]{s), we have 


c(u)f(t + s + u)du 
= jI [0 ,T-t](s) [ e~ {p+q)u {e-KT-t-s-n) _i} du 

b ’ Jo 

= -- I[ 0 ,T-t](s)m(T -t-s ), 

so that 

do(s, f ) = I[o,T-t] («) | - t - s) + /(t + s) | . 



Tims 


(4.6) 


rT—t 


d 2 0 (sj)ds = 


-m{T — t — s) + f(t + s) f ds 


1 

52 


cT-t 


(m(s) + e bs — l} 2 ds. 


For neN and r 6 [0, T — t], we have 


b n (t + u + s, r)f(t + r )dr = -(t> n -\t)e- q ^ a ^m{T - t ). 

0 
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So, the first term on the right-hand side of (14.5p is equal to 

-w^wr (T - t)n 

while the second term to — (1 /b)e~ q ^ t+s ^m(T — t)</> n_1 (t), whence 

2/7P _| Jh+s) 

d n {s, f ) = 4 , 0 M T ~ 


b(p + 2 q) 


Therefore, 


rf ' (S ’ f)dS = b\p + 2qf m2{T ~ l){t) - 


From this and \4>(t)\ < 1, we get 

OO 

E 


“ 2 , , , 2 qe- 2 «‘m 2 (T-t) 

<*„(»>/)<*» = 


(4.7) 


n= 1 


'0 


b 2 (p + 2g) 2 {l — <f 2 (t)} 
2 qm 2 (T — t) 


b 2 {(p + 2 q) 2 e 2qt — p 2 } 

Theorem 13.11 follows from (14. ip , (14.3p , (14.4[) , (14. 6 p , and (14.7p . 


Remark 4.1. If we use other functions as l in (12.21) . then the above proof 
of Theorem, \3.1\ does not work. For, Theorems f.7 and f.12 in [2] are 
no more available. The authors do not know other functions which give 
results similar to those in Section 0 


5. An associated Markov process 


We define 

u(t) := [ e ip+q)s l(s)dW*(s) (t > 0) 

Jo 

so that 

Z(t) = W*(t) - [ pe~ {p+q)s u{s)ds (t > 0), 

Jo 

and couple the short rate process {r(t)} with (w(£)}. Then, the resulting 
two-dimensional process {(r(t), u(t)) T } satisfies the SDE 


r(t) 
u(t) 

. 1 ) 


a 

0 


/b ape / r(t) 

VO 0 ) V u{t) 


dt 


+ ( e<A( t ) ) dW ’V 


(t > 0), 


whence it is a Markov process. This implies that P(t, T), which is defined 
by (12.71) . would have an expression of the form P(t, T) = F(t, r(t),u(t ); T). 
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The aim of this section is to derive an explicit expression for such F and 
to present a numerical method to evaluate contingent claims based on it. 


Lemma 5.1. We have 


a s _ 


P 


,-<i s 


) dZ(s) 


,-pt 


1 - 


pe 


—2 qt 


u(t) (t > 0). 


'o v P+2q~ J ~ ' ' l(t) V" P + 2g, 

Proof. Since w(0) = 0 and e qt dZ{t ) = {e~ pt /l(t)}du(t) — pe~ qt u(t)dt, 


i e -pt rt. r e -ps ^ 1 


u. 


u(s)ds — pe ps u(s)ds. 

Jo 


Similarly, 


p P -(p+2q)t 

/ e~ qs dZ{s) = ——u(t) - 
1 0 L \ l > Jo 


1 ( p-(p+2g)s 'I ' 


l(s) 


u(s)ds 


Combining, 

rt 

,qs 


p 


F qs ) dZ( 


o-pt 


s = 


Jo 
where 


R{s) : = 


p + 2q~ J~ v ~' l(t) 

3 -(P+ 2q)s\ ' 


— pe ( ' p+2qS>s u(s)ds. 

Jo 


pe ? ) u(t)+ f R(s)u{s)ds, 
p + 2qj J 0 


p 


p + 2q 
However, from 


l(s) 


+ pe 


~(p+2q)s 


D -ps 


l(s) 


—ps 


np-ps _ nP -{p+2q)s ( np~(P + 2c d s 

R(s) = - -vf?-+ l ---e- ps 


l(s) 


p + 2q 


l(s) 


pe 


+ p 


and 


np-ip+zq)** 

V — -e" ps 


= 2qe~ ps 


p + 2q ~ J 1 J(s) J [l(s] 

we see that R(s) =0. Thus the lemma follows. 


- 1 


□ 


Here is the explicit expression for F stated above. It may be regarded 
as an affine term-structure for the model (cf. [5], [6] and fTj, Section 5.3]). 


Theorem 5.2. Let P(t,T) be as in Theorem \3.1[ Then we have 
P(t,T) = F(t,r(t),u(ty,T) (0 <t<T), 


F(t, x, y\T ) := exp {-A(t, T ) - C(t, T)x + D(t, T)y} 

((t,x,y) e [0,T] x R x 


where 

(5.2) 
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with deterministic function D(t,T ) defined by 

D(t,T) : = ^e^ p+q)t m(T - t). 

Proof. By Lemma f5.il we have 

U(t, T) = e~ {p+q)t m(T - t)u(t) = D(t, T)u(t). 

This and Theorem 13.11 yield the theorem. □ 


We turn to analogs of term-structure equations for European contingent 
claims in the model AT 

Proposition 5.3. Let S G (0, oo) and let g : M x M —> M and G : [0, S] x 
M x M —* M be continuous. We assume that G is in C 1,2,2 ([0, S) x R x 1) 
and satisfies 

dG 

(5 3) { ~dt^ t,X,V ^ + £ G ( t ’ x >y) = 0 ((T x , y) G [0, S') xlxR), 

G(S, x, y) = g(x, y) ((x, y) G I x R), 


where 


CG := 


2 dx 2 dxdy 


d 2 G {e^H(t)} 2 d 2 G 


dy 2 


dG 


+ {a — bx — pae b+^yX __- X G. 

k J ox 


We also assume 

(5.4) f 

Jo 

where 


Then 


E q 


e -2f* r ( s )ds {AG(t,r(t), U (t))¥ 


ox Oy 


dt < oo, 


(5.5) E q [e-^ r(s)d ^(r(^),u(5)) E t ] = G(t,r{t),u(t )) (0 < t < S ). 

Proof. From (15. 1ft . (I5.3jl and ltd’s formula, we have 

dG(t, r(t),u(t)) = r{t)G{t , r(t),u(t))dt + A G(t, r(t),u(t)))dW*(t ), 
whence 

d{e~ £ r(s)ds G(t, r(t), u(t))) = e~ f ° r(s)ds AG(t, r(t),u(t))dW*(t). 


li 













From this and (15.4p . we see that {e So r ^ ds G(t, r(t),u(t))}o<t<s is a square 
integrable martingale, whence 


e -ti r ^ ds G (t, r (t),u(t)) = E q 

= E q 


~ So r ( s ) ds 


e J o 

'S 


G(S, r(S),u(S)) 


Tt 


e-/o ^) ds g{r{S),u{S )) 


Tt 


Thus (I5.5[) follows. 


□ 


Suppose that, for 0 < S < T and a continuous function h : (0, oo) —> 
R, we want to evaluate a European contingent claim h(P(S,T )) with 
maturity S', written on a T-bond. Since 

h(P(S,T)) = h(F(S,r(S),u(Sy,T)) 

for F in (15.21) . the price of the claim at time t G [0, S] is given by the left- 
hand side of (15.51) with g(x,y) = h(F(S,x,y-,T)). Therefore, by Proposi¬ 
tion 15.31 and under the assumption that (j5.3j) has a suitable solution, the 
evaluation is reduced to solving (15.3|) numerically. In this sense, we may 
regard (15.31) as a term-structure equation for the claim h(P(S,T)) in the 
model M (cf. [H Section 5.2]). 

We illustrate the efficiency of the numerical method stated above with 
two examples. 


Example 5.4. We consider a T-bond. Here we take the values 

T — a = 0.12, b = 1.9, a = 0.35, p = 0.034, q = 0.12. 

In Figure 1, we compare the values of G(0,r(0),0) obtained solving (15. 3 p 
with q(x,y) = 1 numerically and the exact values given by P(0,1) = 
T(0, r(0), 0; 1) with F in (1T2D . 

Example 5.5. We consider a bond call option (. P(S,T ) — K) + with ma¬ 
turity S and exercise price K. We take the values 

T = 1, K = 0.3, a = 0.08, b = 1.5, a = 0.3, 
p = 0.07, q = 0.08, r(0) = 0.025 (or r(0) = 2.5 %). 

In Figure 2, we compare the values of G(0,r(0),0) obtained solving (15. 3 p 
with g(x,y ) = ( F(S,x,y,T ) — K) + and F in (j5.2[) numerically and the 
exact values given by Proposition 13.41 


6. Estimation of parameters 

We estimate the parameters of Ai by fitting the yield curve to actual 
yield data. More precisely, we search for the values of parameters such 
that the following least squares error is minimized: 

n 

J2{y(0.T,)-Y(0,T t )} 2 

i= 1 
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(see [3 Chapter 3] and P Section 30.8] for related discussions on cali¬ 
brating term-structure models). Here, Y(0,T) is the yield for a T-bond 
in the model Ml given by (13.11) . and t/(0,Tj)’s are actual yield data. We 
take n = 10 and 

13 6 

T\ = —, To = —, T 3 = —, T 4 = 1 , T 5 = 2 , 

1 12 12 12 

T 6 = 3, T 7 = 5, T 8 = 7, T 9 = 10, T 10 = 20 (in years). 

As 2 /( 0 , Tj)’s, we use the Daily Treasury Yield Curve Rates published at 
the US Treasury Website. 

First, we use the data on December 31, 2007. The estimated values of 
the parameters are as follows: 

a = 0.1635, b = 1.8952, a = 0.7247, p = 0.0909, q = 0.2100, 

r(0) = 0.0240 (or r(0) = 2.40 %). 

In Figure 3, we show the data y( 0, T*) and the fitted yield curve Y (0, T). 
For comparison, we also show the Vasicek yield curve fitted to the data. 

Next, we use the data on May 24, 2005. The estimated values of the 
parameters are as follows: 

a = 0.0822, b = 1.5561, a = 0.3007, p = 0.0696, q = 0.0758, 

r(0) = 0.0259 (or r(0) = 2.59 %). 

The data and fitted yield curves are shown in Figure 4. 

Finally, we use the data on September 15, 2008. The estimated values 
of the parameters are as follows: 

a = 0.1216, b = 1.6806, a = 0.6246, p = 0.1170, q = 0.1623, 

r( 0 ) = 1.0010 x 10“ 5 (or r( 0 ) = 1.0010 x HT 3 %). 

The data and fitted yield curves are shown in Figure 5. 
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Figure 3. Yield Curve on December 31, 2007 
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15 






4.5 



Figure 5. Yield Curve on September 15, 2008 
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